Competition— induced shifts of spectral peaks in photon coincidence spectroscopy 
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We show that shifts in locations of two-photon coincidence spectral peaks, for a bichromatically- 
driven two-level atom passing through a single-mode cavity, are due to competition between excita- 
tion pathways for a Jaynes-Cummings system. We also discuss an analogous shift of (single-photon) 
spectral peaks for a driven three-level V-system, which demonstrates that competition between ex- 
citation pathways is also important in this simple system. 



PACS numbers: 42.50.Ct, 42.50.Dv 
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I. INTRODUCTION 

Cavity quantum electrodynamics (CQED) has contin- 
ued to develop rapidly, driven both by recent experimen- 
tal successes and by the promise of exciting new appli- 
cations. Advances in atom cooling techniques, as well as 
development of high-Q optical cavities with large-dipole 
coupling, have enabled testing of the strong-coupling 
regime of CQED Q as well as the trapping of single 
atoms in optical cavities 0, Q . 

Applications of CQED, especially for such applications 
as the quantum logic gate jj| , rely critically on entangle- 
ment between the field degree of freedom and the inter- 
nal electronic state of the atom |^ 0| . This entanglement 
is not only challenging to achieve, it is also difficult to 
probe. In the optical regime of CQED, photon coinci- 
dence spectroscopy (PCS) has been proposed as a feasible 
and unambiguous method for detecting genuine quantum 
phenomena in CQED. This technique employs a bichro- 
matic (or multichromatic) driving field acting on the 
combined atom-cavity system and detects two-photon 
(or multiphoton) decays, respectively, by registering pho- 
ton coincidences in the cavity field emission (H Q, H • 

The simplest case of PCS is two-photon coincidence 
spectroscopy (2PCS) for probing the nonlinear portion 
of the Jaynes-Cummings (JC) spectrum. The technique 
of 2PCS proceeds, first by driving the atomic beam with 
a bichromatic field (consisting of a fixed driving field with 
fixed frequency wi and a scanning field with tunable fre- 
quency L02), which causes two-photon excitation to the 
second couplet of the JC ladder, followed by two-photon 
decay from the atom-cavity system (Fig.^. The objec- 
tive is to count photon pairs emitted from the cavity as 
a function of varying 0^2 • When the sum frequency for 
the bichromatic driving field oji +L02 matches a transition 
frequency from the second couplet in the JC spectrum to 
the ground state, the result is an enhanced two-photon 
count rate (2PCR). 

This simple picture of excitation pathways agrees with 
simulations [a, Q , but displacement of coincidence peaks 
as a function of scanning field frequency was evident yet 
not properly understood Q\ . Here we show that the devi- 
ation of peak shifts is due to competition between excita- 
tion pathways. This picture yields excellent quantitative 



agreement with simulations of 2PCS. 



II. MASTER EQUATION AND TWO-PHOTON 
COUNT RATES 

In the electric-dipole and rotating-wave approxima- 
tions, the JC Hamiltonian for the two-level atom (2LA) 
coupled with the single mode is 



H{g) = Lo{(Jz + o)a) + ig{r){a) 



acT+} 



(1) 



with r the position of the atom, g{r) the position- 
dependent dipole coupling strength, a and the anni- 
hilation and creation operators for photons in the cav- 
ity field, <J-|-, (T_, and ct^ the 2LA raising, lowering and 
inversion operators, respectively, and h = \. Provided 
that the atoms move sufficiently slowly through the cav- 
ity IE Q 1 the atom can be treated as if it were at rest 
at some randomly located position r. As the position r 
is a randomly varying quantity, the value of the cou- 
pling strength g itself is also random. Hence, a cou- 
pling strength distribution P{g) can be constructed l7|, 
and we assume the P{g) depicted in Fig. 5 of Ref. 0. 
This P[g) corresponds to the coupling strength distribu- 
tion obtained for atoms passing through a rectangular 
mask, centered at an antinode, for the cavity sustaining 
only the TEMqo mode. The dimension of the mask is wq 
(cavity mode waist) by A/10 (with A the optical wave- 
length). We restrict F^max < 9 < 5max for gmax the 
coupling strength at an antinode along the cavity longi- 
tudinal axis and F an effective cut-off term. 

The spectrum for the Hamiltonian JQ), depicted in 
Fig. n is the well-known JC spectrum, or 'ladder' 1^. 
The 'dressed states' of the combined atom-cavity sys- 
tem are here designated by the lowest-energy state |0) = 
|0)cav <8) |g)atom = |0,g), and, for n a positive integer, 
the 'excited' states \n)± = i/\/2 (|n — 1, e) ± z |n, g)), 
with \n) the Fock state of the cavity mode and |g) (|e)) 
the ground (excited) state of the 2LA. The effect of aver- 
aging over P{g) is an inhomogeneous spectral broaden- 
ing, due to atomic position variability, as shown in Fig. ^ 

Two-photon excitation is provided by driving the atom 
directly, as it traverses the cavity, with a bichromatic 
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FIG. 1: Two-photon excitation scheme from the ground 
state |0) to the first two excited couplets |n)^ (n £ {1,2}, 
e £ { — ,+}) of the dressed states. The inhomogeneous broad- 
ening of the couplets 1 1) and |2) is 2hgma,x and 2\^hgma.x, re- 
spectively. Two two-photon excitations to the second couplet 
are depicted for a bichromatic driving field with one compo- 
nent of amplitude £i and the other with amplitude £2- The 
excitation pathway on the right {iOi then UJ2) excites reso- 
nantly from |0) to and then may excite resonantly to 
either |2)±. The excitation pathway on the left {uj2 then uji) 
excites resonantly from |0) to to |2)_ for g = {\/2— l)gf- 



system from the ground state |0) to the excited state 

for the subensemble g = g f = u — uji; this subensemble 

corresponds to P{g) = 5{g — gj). 

The scanning-field frequency is lj2^ and excites the 
subensemble for g = gf from to one of the two 

states in the second couplet of the JC ladder, namely 
|2)±. Thus, the range of scanning frequencies for uj2 
must include the |1)_ < — > |2)-|- transition frequencies, 
LJ ± (\/2 T 1)3, respectively. The amplitudes of the two 
chromatic components should be large enough to ensure 
sufficient occupation of the excited state but not large 
enough that significant Stark shifting or nonnegligible 
occupation of the higher-order states occurs. Enhanced 
rates of photon pair detection are then sought as the 
scanning frequency UJ2 is varied such that uji + UJ2 is res- 
onant with some transition |0) < — > |2)± as depicted in 

Fig.m 

The Born-Markov approximation is applied to both 
radiation reservoirs: the reservoir for the field leaving the 



cavity and the reservoir for direct fluorescence of the 2LA 
from the sides of the cavity into free space. The cavity 
damping rate is k, and the cavity-inhibited spontaneous 
emission rate into free space is 7. The master equation 
for this system Q can be expressed as p = Cp for C the 
Liouville operator. Here C = £cff + f + i/, i.e. a sum 
of a Liouville operator >Ccff, £^ri explicit time-dependent 
Liouville operator T) and a 'jump' term J7. We introduce 
5 = LO2 — LO1 and work in the rotating picture with respect 
to the driving-field component wi. 

The effective Hamiltonian is (without jump terms) 

H,tf{g,£i) = {oj-LUi){a,+a^a) + Eig) + T{£i) 

— iKfl^a — i(7/2)(7_|_(T_, (2) 



for 2(17) — ig{a^a- — aa 

ator and T{£i) — i£i{a^ — ct_) a monocuromaiic zi 
driving term. The corresponding Liouville operator is 



the quantum exchange oper- 
^ monochromatic 2LA 



^cs{g,£i)p = Hcff{g,£i)p- pHlff{g,£i 



(3) 



The operator for the jump term is such that 

J^p ~ 2fiapa^ + ja^pa^, (4) 

and the time-dependent Liouville operator is 'Dit) such 
that 



Vp=-t [T{£2e-'''),p] 



(5) 



Solving the master equation for p, and averaging over 
P{g), which accounts for atomic position variability, 
yields a solution 



P{9)p{9)dg. 



(6) 



For a bichromatic driving field, the density matrix does 
not settle to a steady state value, but in the long-time 
limit t — > 00, the Bloch function expansion is 



lim p{t) 

r— i-OO 



m— — 00 



ArnSt 



(7) 



with p„j time-independent matrices. As the photocount 
integration time is expected to be long compared to the 
detuning S, it is reasonable to assume that rapidly oscil- 
lating terms average out and therefore approximate p{g) 
by truncating the expansion Q. 

The optical signature for entanglement is the 2PCR 



w 



(2) 



{S,£,) = {:n' ■.){S,£,), 



(8) 



for n = a^a. The 2PCR is shown in Fig.[21as a dotted line, 
and the peaks expected a.t 6 = {6 — gf)/gf — ±(V2 — 1) 
are not distinguishable. This indistinguishability is due 
to off-resonant two-photon transitions with both pho- 
tons of frequency W2, and may be overcome via the 
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FIG. 2; The 2PCR vs the normalized detuning S for the 
inhomogeneously broadened system, with £i/k. — l/\/2, 
£2/1^ = \/2, = 9, 7/k — 2: to'-^'((5, £1) is the soUd hne 

and A^^\5,£i) is the dotted Une. 

method of background subtraction: the experiment is re- 
peated with £1=0 and the resuhant 2PCR subtracted. 
Thus, we define a diff'erence-2PCR 

A'^^\S,£,) = w^^\S,£i)-w^^HS,£i = 0). (9) 

For a specific set of parameters we show in Fig. El the 
2PCR with and without background subtraction, which 
is significant in the domain of homogeneous broaden- 
ing — 1 < ^ < 1, but the region near 5 = 1-1- is not 
significantly affected by background subtraction. Back- 
ground subtraction helps to resolve the two two-photon 
spectral peaks at 6 — ±(-\/2 — 1). 

III. COMPETITION-INDUCED SPECTRAL 
PEAK SHIFTS 

The objective of 2PCS is to obtain a 2PCR with peaks 
at values of S which cannot be explained by semiclassical 
models. The peaks at 6 = ±(-\/2 — 1) and S — \/2 -t- 1 
involve the quantity \/2, which follows from the quantum 
field expression a|2) = V2|l)- This V2 is therefore a gen- 
uine, unambiguous quantum field signature. However, 
the actual peaks in 2PCR cannot be expected to occur 
at exactly these values of 6 as observed quite clearly in 
Fig. 5 of Ref. In Fig. Owe show in detail the peak 
near S = \/2+l as a solid line for P{g) — (>{g — gf) and an 
anomalous shift from its expected position at S = \/2 + l, 
or A5 = S~ (\/2+ 1) = 0, is evident. We begin exploring 
this anomalous shift for the g — gf subensemble. 

We observe that this peak is noticeably shifted from 
the expected location. The origin of the shifts in peaks is 
subtle. In Ref. 01 the suggestion is that "the peaks are 
shifted slightly, probably due to the motion of atoms". 



Although atomic motion may contribute to peak shifts, 
we show in this section that competition between exci- 
tation pathways is the major factor in shifts of 2PCS 
peaks. 

The peak at 5 = \f2 + 1 is the most interesting of the 
three major 2PCR peaks (i) because background sub- 
traction is not required to detect this peak (as 8 lies out- 
side the region of inhomogeneous broadening) and (ii) 
because the peak is quite distinct from the large neigh- 
boring pair of peaks at 6 — ±(-\/2 — 1). The excitation 
pathway to obtain the peak at ^ = \/2 + 1 is via an l>j\ 
photon which resonantly excites |0) < — > followed 
by an L02 photon which resonantly excites |1)_ < — s- |2)-)_ 
for the 2PCR peak at b = \/2 -f- 1 as shown in Fig. ^ 
However, level |2)_ becomes populated as well as |2) + 
by off-resonant excitation. 

The competition between populating levels |2)± is ev- 
ident in Figs. 3(d,e) of Ref. which depict the popu- 
lations of levels |2)± as functions of g and <5. The |2)_ 
level is noticeably occupied for all values of b for g in the 
vicinity of gy. This occupation of |2)_ was not explicitly 
discussed in Ref. 0, Q , and the significance of this popu- 
lation of |2)_ was not provided with respect to the shift 
in the peak at 8 — y/2 + 1. The shift in the peak is shown 
but not discussed in Ref. . The following analysis pro- 
vides the first detailed picture for explaining the shift in 
the 2PCR peaks. 

The effect of competition is made clear in Fig.O Here 
we present one portion of the 2PCR peak as a function 
of 5, in the vicinity of AS = 0. The solid line corresponds 
to the 2PCS. The dotted line corresponds to the simu- 
lated 2PCS but with the artificial restriction that the 
matrix elements of the Hamiltonian driving field opera- 
tor T(£2 exp(— i(5i)) for the transition |1)_ < — > |2)_ are 
set to zero. That is, we artificially impose the condition 

_(l|T(£2expH,5i))|2>- = =_(2|T(£2expH,5i))|l>- • 

(10) 

By doing so we observe, and show in Fig.|3| that the mod- 
ified 2PCS graph is indeed centered at A5 = (equiva- 
lently, at 1=^/2 + I). The comparison between the two 
cases, with both _(1|T(£2 exp(— i(5t))|2)_ and its con- 
jugate set to zero vs both not being set to zero, pro- 
vides convincing evidence that the shift in the peak is 
indeed primarily due to competition between two excita- 
tion pathways. 

Thus far our attention has been restricted to the 
P{g) — S{g~ gf) case. This case is a significant contribu- 
tor to the peak shift for the inhomogeneously-broadencd 
peak shift obtained in the full simulation, accounting for 
atomic position variability. The simulation of the 2PCR 
peak near 6 = -I- 1 is depicted in Fig. 0{a) as a solid 
line. This solid line evidently has a peak at AS — —0.141. 
The above reasoning, for the P{g) = S{g — gj) case, sug- 
gests that the competing pathway |1)_ < — > |2)_ is par- 
tially responsible for this peak shift. Therefore, we re- 
peat the simulation, but with the restriction (|10|l . The 
result of this simulation is also depicted in Fig. EJa) as 
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FIG. 3: Two-photon count rate (2PCR) vs the shift in nor- 
malized scanning frequency of the peak at S — \f2 + 1 
for gjn = 9, fi/fi; = 1/^/2, £2 Ik. = y/2 and 7/^ = 2 de- 
picted as the solid line. The dotted line corresponds to the 
case with condition lllUl imposed. 



a dashed line. The new peak position is located ap- 
proximately at A6 = —0.094. This peak position can 
be shifted slightly by eliminating other pathways. By 
fixing all transition terms of the type H10|l to be zero 
except for (0|T(f i)|l)_, _ (1|T(£2 exp(-i(5i))|2) + , and 
their complex conjugates, the resultant peak is centered 
at AS = —0.091. These additional pathways, which were 
ignored in obtaining A6 = —0.094, do not contribute 
much to the overall peak shift, and we do not depict this 
specific peak shift. 

To account for the remainder of the peak shift observed 
in Fig.^a), we determine other pathways to an enhanced 
2PCR, which arise in the full inhomogeneously broadened 
case but not in the g ~ gf case. We now consider the ef- 
fect of the homogeneous linewidth at the level state. 
To incorporate this effect, we solve the dynamics for the 
effective Hamiltonian (j^J plus the driving term ^ but do 
not include the jump term I^J. The 2PCR, which elim- 
inates the pathway < — > |2)„ and the homogeneous 
linewidth at is depicted in Fig. Efb). We observe 
that the resultant peak shift is now A5 = —0.057. The 
consideration of the homogeneous linewidth does matter 
in accounting for the peak shift. 

To understand the peak shift due to the homogeneous 
linewidth at we consider the following. The path- 

way |0) < — > < — > |2)+ is relatively unaffected by 
the homogeneous linewidth as it is a resonant transition 
at AS = 0, but the homogeneous linewidth does enable 
off-resonant excitation to |2)_, even after the restriction 
()10|l has been applied. The homogeneous linewidth of 
the |1)_ level enables two-photons transitions to take 
place for a range of coupling strength g, via the sequence 
of an uji photon and then an uj2 photon. By eliminat- 
ing this homogeneous linewidth, the loi photons really 
must satisfy the condition lu ~ uji = g f in order to off- 
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FIG. 4: The 2PCR vs peak shift from 5 = \/2 -f 1 for 
the inhomogeneously broadened system with £i/k = l/\/2, 
£2/^^ = V2, Qf/K = 9, 7/« = 2. Two cases are shown: 
(a) with (solid) and without (dots) the imposition of con- 
dition Hl()|l : and (b) with the imposition of condition HIO^ and 
the homogeneous linewidth (and jump terms) set to zero for 
the 1 1) _ level of the JC ladder (solid) , and with the additional 
imposition that {0|T(f2 exp(— i(5t))| 1)+ and its complex con- 
jugate are set to zero (dots). 
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resonantly excite to the |2)_ level. The homogeneous 
linewidth is responsible for permitting competing path- 
ways from |0) < — > < — > |2)_ for a range of g values, 
and eliminating this homogeneous linewidth eliminates 
a whole class of competing pathways for the inhomogc- 
neously broadened (general P{g)) system. 

The final significant contribution to the peak shift of 
the 2PCR curve arises for small g. The low coupling 
strength limit is important because the distribution P{g) 
is so heavily weighted in favor of low g. When the cou- 
pling strength is small, the splitting between the first 
couplet states |1)± becomes correspondingly small. Con- 
sequently one would expect a competing pathway to arise 
via excitation to + . We eliminate the contribution due 
to the pathway |0) < — > Incorporating this restric- 

tion with the other restrictions previously mentioned 
(eliminating the < — > |2)_ pathway and the homo- 
geneous linewidth at leads to the dotted curve in 
Fig.lHb). The resultant peak shift is just AS -0.018. 

The three effects considered above account for 87% of 
the peak shift. Better agreement can be obtained by re- 
ducing homogeneous linewidth on the |2)_ level or by 
eliminating transitions to the competing state |2)_ alto- 
gether. However, the agreement demonstrated here ex- 
hibits excellent agreement with the full simulated 2PCR 
peak at 5 = \/2 + 1. There is no doubt that this shift is 
primarily due to competing pathways in the JC ladder. 

In accounting for the peak shift due to competing path- 
ways, we have included restrictions of coherent transi- 
tions and thereby accounted for the majority of the peak 
shift. There is also competition with incoherent path- 
ways, which are accounted for by the quantum jump 
terms. The quantum jumps correspond to the transfer 
of quanta between levels of the JC ladder, but the con- 
tribution of these jumps to peak shifts is not large for 
the bichromatically-driven JC system. Never the less, 
the jump terms are responsible for part of the peak shift. 
This competition with incoherent pathways is quite pro- 
nounced for the monochromatically-driven three-level 
system discussed in the next section. 



IV. THE MONOCHROMATICALLY-DRIVEN 
THREE-LEVEL V-SYSTEM 



troscopy. The normalized frequency of this driving field 
is (5 = ~^/9f for S = uj — LJi. We expect a peak at 6 — 1, 
for which resonant excitation occurs if we ignore compe- 
tition between excitation pathways. 

In Fig.|5lwe see that, for A6 — S—1, the peak is shifted 
to the left. Competition due to off-resonant excitation 
is established by setting certain matrix elements of the 
driving term in the master equation to zero, by analogy 
with the 2PCR considered above. That is, for Hcs{g, £i) 
in Eq. (|2Jl, we impose 

(0|T(£:i exp(-i<5t))|l)_ = =-{l\T{£iexp{-iSt))\0) 

(11) 

in the simulation. We observe that this condition pro- 
duces a plot of {a^a) that is exactly centered at AS = 0. 
Therefore, the competition with the off-resonant excita- 
tion pathway |0) < — > |1)_ is entirely responsible for the 
peak shift of (a^a) vs S. 

Another perspective for understanding the spectral 
peak shift for the driven V-system is obtained directly 
from the effective Hamiltonian of Eq. Q solved in the lin- 
ear regime. We define for the pseudo-pure state V', Cq = 
(OjV') and Ci± =± (lIV")- The resulting set of differen- 
tial equations for these dressed-state amplitudes Cq, Ci_ 
and Ci+ is given by 



Co = 

Ci- = 

Ci+ = 



^(Ci_ + Ci+), 



(12) 



1 



i)Cr. 



(13) 



*('5 + 5) + ^(«+|) 
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(14) 



Taking the Laplace transform in the long time limit, a 
minimum and two maxima occur in the PCR for 



Sq — (minimum) 



(15) 



/ 2 -I- £2 

S = d± = —\lzL (maxima), (16) 

V ^9 



The shift in spectral peaks should not be unique to 
a bichromatically-driven JC system. A simple case to 
consider is the monochromatically-driven three-level V- 
system (the JC ladder in the linear regime, cf Fig. 
which is of interest for studying normal-mode (or vacuum 
Rabi) splitting [TSIHIil. 

We consider the three-level truncation of the JC ladder 
to {|0), |1)±} and a monochromatic driving field which is 
resonant with the |0) < — > |1)+ transition {uji = uj + gf) 
where tui is allowed to vary. In contrast to the previ- 
ous section, enhanced one-photon count rates identify 
resonances and correspond to the standard method of 
spectroscopy, as opposed to photon coincidence spec- 



for the decay rate j/k = 2. To identify competition be- 
tween pathways, we impose the condition Eq. Hll|) . This 
condition modifies Eq. H12|l by replacing the Ci_ term 
by zero and modifies Eq. 1)1 3|l by replacing Cq by zero. 
The peak at (5_ thus vanishes, as there is no driving to 
the |1)_ state. The peak at d+ is shifted to ^ = 1 (where 
the peak would be expected initially without any com- 
petition between excitation pathways). The formula for 
the competition-induced peak shift for the three-level JC 
system is thus given by 



AS = 6,-1 = X 1- 



2 + g2 



(17) 
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FIG. 5: One-photon count rate (IPCR) of a three-level sys- 
tem vs the shift in normalized scanning frequency of the 
peak at (5 = 1 for g/n = 9, Si/k, = \/2 and 7/k; = 2 de- 
picted as the solid line. The dotted line corresponds to the 
case {Q\T{£iexp{-i5t))\l)- = (l|T(£:i exp(-i(5t))|0). 
The vertical line corresponds to A J in Eq. I|17|l . 



The peak shift AS presented in Eq. H17|l accounts par- 
tially for the total peak shift presented in Fig. The 
discrepancy is due solely to the effect of quantum jump 
terms, which have been ignored in Eqs. H15|l - ()16|l . The 
necessity of including quantum jump terms arises be- 
cause dissipation terms cannot be dissociated from fluc- 
tuations in a proper analysis. The quantum jump terms 
are responsible for incoherent transitions of quanta from 
one level to another of the three-level system. We ob- 
serve, in Eqs. ((T^ - fTHjl . that incoherent transitions are 
not included in the mathematics. Specifically, we can 
see that diagonal matrix elements poo = C^Cq and 
pf^ = C*_|_Ci± are not coupled together. This cou- 
pling of diagonal elements of the density matrix arises 
via quantum jump terms. These quantum jump terms 
are responsible for the discrepancy between the pre- 
dicted peak shift given by Eq. ifTTIl and the observed 
peak center of the solid locus in Fig. [3 Whereas, for 
the bichromatically-driven JC system, the contribution 
of incoherent pathways to the net peak shift is small, 
the contribution of incoherent pathways to the peak 
shift for the monochromatically-driven three-level V sys- 
tem is relatively large because there are fewer coherent 
pathways to producing a peak shift as compared to the 
bichromatically-driven JC system considered in the pre- 
vious section. 

The advantage of using Eqs. (|15|I - H16|) . to confirm that 
the peak shift is indeed partially due to competition be- 
tween coherent excitation pathways, is that this method 
is entirely analytical. The analytical method provides 
a simpler conceptual framework to understand the role 
of both coherent and incoherent excitation pathways. 
In this analytical method, we obtain solutions via the 
Laplace transform method rather than relying on nu- 



merical solutions to the full master equation. Numeri- 
cal methods are necessary, though, to fully account for 
spectral peak shifts. 



V. CONCLUSION 

The analysis presented here provides convincing evi- 
dence that the shifts in two-photon coincidence spectral 
peaks, from expected locations suggested by the JC lad- 
der of Fig. ^ are due overwhelmingly to competition be- 
tween excitation pathways. 

We have analyzed in detail the two-photon spectral 
peaks, outside the domain where background subtrac- 
tion is essential for resolution. The shift of the spectral 
peak has been shown to be a consequence of competition 
between excitation pathways. This demonstration relies 
on numerical methods whereby the competing pathway 
have been artificially eliminated by imposing the condi- 
tion that certain matrix elements are zero in the master 
equation. 

The key contribution to the peak shift were identi- 
fied to be competition with the |0) < — > < — > |2)_ 
pathway for g « g/, |0) < — > < — > |2)± for 5 ~> 
and the homogeneous linewidth about which al- 

lowed off-resonant excitation to |2)±. Eliminating these 
three effects reduces the two-photon coincidence peak 
shift by 87%. 

These shifts in spectral peaks are not restricted to 
2PCS. To verify this, we considered the shifts of spec- 
tral peaks for standard (one-photon) spectroscopy and a 
three-level monochromatically-driven V-system. Such a 
system corresponds to the driven JC system in the linear 
regime, which is of interest in the study of normal-mode 
(vacuum Rabi splitting). We also observe the shift of 
spectral peaks in our simulations. By ignoring the jump 
terms, the equations become easily solvable and we verify 
that a significant part of the peak shift for the V-system 
are indeed due to competition between excitation path- 
ways. The analytical method enables us to clearly see the 
distinction between coherent and incoherent pathways as 
competition-induced peak shifts. The contribution of 
quantum jump terms to the net peak shift is significant 
for the monochromatically-driven V system, whereas the 
effect is much smaller for the bichromatically-driven JC 
system. The combination of competition with coherent 
and incoherent pathways fully accounts for the peak shift 
of the monochromatically-driven three-level V system 
analyzed in Sec. II VI 

These results clarify the expected peak locations in 
2PCS and should also apply to the more general Mul- 
tiphoton Coincidence Spectroscopy 0. As photon coin- 
cidence spectroscopy must be accurate to verify the sig- 
nature of entanglement, these peak shifts must be fully 
accounted for and properly understood. The analysis 
presented here does indeed explain the peak shifts clearly 
and quantitatively. 
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